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Abstract. We highlight modular equations discovered by Somos and Ramanu- 
jan, and use them to prove new relations between lattice sums and hypergeometric 
functions. We also discuss progress towards solving Boyd's Mahler measure con- 
jectures, and conjecture a new formula for L{E, 2) of conductor 17 elliptic curves. 



1. Introduction 

Modular equations appear in a variety of number-theoretic contexts. Their con- 
nection to formulas for I/tt [14], Ramanujan constants such as e'^^'^^ [21], and 
elliptic curve cryptography is well established. In the classical theory of modular 
forms, an nth degree modular equation is an algebraic relation between j(r) and 
j(nr), where j(r) is the j'-invariant. For our purposes a modular equation is simply 
a non-trivial algebraic relation between theta (or eta) functions. In this paper we 
use modular equations to study four-dimensional lattice sums. The lattice sums are 
interesting because they arise in the study of Mahler measures of elliptic curves. 

There are a large number of hypothetical relations between special values of L- 
series of elliptic curves, and Mahler measures of two-variable polynomials. The 
Mahler measures m{a), n{a), and g{a) are defined by 



ma 



[ \og\y + y^^ + z + z-^ + a\d9id92, 
Jo 

(1) n{a) := / \og\y^ + + 1 - ayz\ d9id92, 

Jo Jo 

g{a):=[ [ log \{y + 1) {z + l){y + z)- ay z\d9,d92, 
Jo Jo 

where y = e^'^*^^, and z = e^'^*^^. Boyd conjectured that for all integral values of 
7^4 [6]: 

m{k) = \LiE,2), 

where E is an elliptic curve, q is rational, and both E and q depend on k. He also 
discovered large numbers of formulas involving g{a) and n{a). In cases where E 
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has a small conductor, it is frequently possible to express L{E, 2) in terms of four- 
dimensional lattice sums. Thus many of Boyd's identities can be regarded as series 
acceleration formulas. The main goal of this paper is to prove new formulas for 
the lattice sum F{b, c), defined in (12). So far there are at least 18 instances where 
F{b, c) is known (or conjectured) to reduce to integrals of elementary functions. The 
modular equations of Somos and Ramanujan are the main tools in our analysis. 

2. Eta function product identities 

Somos discovered thousands of new modular equations by searching for linear 
relations between products of Dedekind eta functions. Somos refers to these formulas 
as "eta function product identities" . The existence of eta function product identities 
can be established by using fact that certain modular parameters (such as j(r)) equal 
rational expressions involving eta functions. By clearing denominators it is possible 
to rewrite classical modular equations as eta function product identities. A major 
surprise of Somos's experimental approach, is that it turned up a large number of 
unexpectedly simple identities. In order to give an example, first consider the eta 
function with respect to q: 

oo oo 

V{q) = q'"' 11(1 -q-)=Y. {-iTq^'-^'^' 

n=l n=—oo 

and adopt the short hand notation 

The following formula is the smallest eta function product identity in Somos's list 
[19]: 

(2) 6266610630 = 61612615620 + 636465660- 

Notice that all three monomials are products of four eta functions, and are essentially 
weight-two modular forms. There are no known identities between eta products of 
weight less than two, and (2) appears to be the only three-term linear relation 
between products of four eta functions. Many additional identities are known if 
the number of terms is allowed to increase, or if eta products of higher weight are 
considered. For additional examples see formulas (14), (22), (23), (26), and (29). 

Identities such as (2) can be proved almost effortlessly with the theory of modular 
forms. A typical proof involves checking that the first few Fourier coefficients of 
a presumed identity vanish. Sturm's Theorem furnishes an upper bound on the 
number of coefficients that need to be examined [13]. We note that it is often 
possible, but usually more difficult, to prove such identities via elementary g-series 
methods. Ramanujan was a master g-series manipulator, and his notebooks are 
filled with various modular equations and their corollaries. We conclude this section 
by providing a Ramanujan-style proof of (2). The main news is that (2) can be 
derived from modular equations known to Ramanujan. 

Theorem 1. The following formula is true: 

(3) 6266610630 = 61612615620 + 636465660- 
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Proof. Before proving (3) we need to define a small amount of notation. Let us 
denote the usual theta functions by 

oo oo 

(4) ^{q):=Y,q-\ ^'(g) := E 

n=— oo n=0 

Furthermore define Uj and zj by 

Uj . i ^^^^.^ , . if [q ). 

Notice that Ramanujan uses a slightly different notation [2]. He typically sets a = 
Ui, and he says that "/3 has degree j over a" when /3 = Uj. Certain values of the 
eta function can be expressed in terms of Ui and zi [2, p. 124]. We have 

(5) r^(g)=2-Vy/^^(l-«0^/^v^, 

(6) v{q') =2-'/-' {u,{l - u,)}'/'' 

(7) v{q')=2~'/'u/\l-u,y/''^,. 

Now we prove (3). By (6) the left-hand side of the identity becomes 

6266610630 = 2~'^/^{miM3M5Mi5(1 - Mi)(l - M3)(l - U^)^! - Uir^)}^/^"^ Z1Z3Z5Z15. 

By (5) and (7), the right-hand side of the identity becomes 
61612615620 + 636465660 

=2-5/3 ({^3^5(1 - ni)(l - ni5)}^/^ {nini5(l - n3)(l - u,)}'^'^ 

+ {MiMi5(1 - M3)(l - ^5)}^^^ {^43^5(1 - Ui){l - ^15)}^^^'^) y/z^Z^Z^. 

Combining the last two formulas shows that (3) is equivalent to 

2^/'^ {^1^3^5^15(1 - ni)(l - n3)(l - n5)(l - ^15)}^/^^ 

= {m3M5(1 - Mi)(l - ^15)}^^^ + {miMi5(1 - - ^^5)}^^*^ • 

It is sufficient to show that (8) can be deduced from Ramanujan's modular equations. 

The first modular equation we require can be recovered by multiplying entries 
11.1 and 11.2 in [2, p. 383]: 

((«i«i5)'/' + {(1 - «i)(l - «i5)}'/') [{u.u.y/' + {(1 - «3)(1 - u,)}'/') 
Rearranging yields an identity for the right-hand side of (8): 

{U3U5(1 - - ^15)}'/' + {^^1^^15(1 - ^3)(1 - ^5)}'/' 

= 1 - {1x1^3^5^15}^/' - {(1 - ni)(l - U3)(l - n5)(l - ni5)}'/^ 
By entry 11.14 in [2, p. 385], it is clear that 

1 - {^1^31151115}^/" - {(1 - ni)(l - n3)(l - n5)(l - ni5)}^/' 

= 2l/3{nili3li5lii5(l - lii)(l - li3)(l - 115)(1 - lii5)}^/'^ 



1. 
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The theorem follows from combining (9) and (10) to recover (8). □ 

It appears to be very difficult to explain why identities such as (3) exist. Our ini- 
tial motivation for constructing an elementary proof of (3), was to find a method for 
generating more identities. It would be interesting if a systematic method for gen- 
erating weight-two eta product identities could be discovered. This is a surprisingly 
important question for studying lattice sums and the Beilinson conjectures. 

3. Lattice Sums 

In this section we investigate four-dimensional lattice sums. Many of these sums 
appear in the study of Mahler measures of elliptic curves. Let us define 

F(a, 6, c, d) :=(a + b + c + df 

^ ]^'^ni+n2+n3+n.4 



X 



(a(6ni + 1)2 + 6(6^2 + 1)^ + c(6n3 + 1)^ + diQm + 1)2)" 

The four- dimensional series is not absolutely convergent, so it is necessary to employ 
summation by cubes [5]. Notice that Euler's pentagonal number theorem can be 
used to represent -F(a, 6, c, d) as a useful integral 

(11) F(a,6,c,rf) = J<^ + b + c + df /'\(^a)^(^^.)^(^c)^(g^)iogg^. 

We also use the shorthand notation 

(12) F{b,c) :=F(l,6,c,fec), 

since we are primarily interested in cases where a = 1, d = be, and b and c are 
rational. 

The interplay between values of F{b,c), Boyd's Mahler measure conjectures, and 
the Beilinson conjectures is outlined in [16]. If (6, c) G and (1-1-6) (l-|-c) divides 24, 
then F{b, c) = L{E, 2) for an elliptic curve E. Formulas are now rigorously proved 
relating each of those eight cases to Mahler measures such as m{a) [22]. Because 
Mahler measures often reduce to generalized hypergeometric functions, many of 
Boyd's identities can be regarded as series transformations [15], [12]. It is known 
that 



m{a) = Re 
n(a) = Re 



2 

log(a) -4F3 



9 



2 

log(a) ^4^3 



3 





16 


2,2,2 


a2 


f ^ ^11 


27 


2,2,2 








1/3 i ' 


if |( 



, if a 7^ 0, 

, if |a| is sufficiently large. 



3g{a) =n y — ^/Tj ~'~ y 1/3 J ' 1*^1 sufficiently large. 

The function m{a) also reduces to a 3F2 function if a G M [11], [16]. The first author 
and Zudilin [17] recently proved that 

47r2 7r2 / 1 1 1 1 

(13) F(3,5) = m(l) = — 3F9 ''l^'' — 

^ ' ^ ' ^ 15 ^ ^ 15 V I'i 16 
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Equation (13) is equivalent to a formula that was conjectured by Deninger [8], and 
which helped motivate Boyd's seminal paper [6]. One of the main results in [16], 
is that it is also possible to find formulas for values such as F(l,4) and F{2,2). 
These cases are probably not related to elliptic curve L-values. As a result it was 
hypothesized that it should be possible to "sum up" F{b, c) for arbitrary values of 
b and c. 

3.1. Lacunary cases. It is typically very difficult to prove formulas such as (13). 
The proof of (13) is a g-series proof which utilizes the integral representation (11). In 
general the difficulty of dealing with a lattice sum depends on whether it is lacunary 
or non-lacunary. The lacunary values can be reduced to two-dimensional sums, 
which are (almost) always easier to deal with than four- dimensional sums. It is 
often difficult to determine whether or not a particular sum is lacunary. Cases such 
as F(l, 1), F(l, 2), and F(l, 3) can be equated to L-values of CM elliptic curves, and 
their lacunarity follows from the CM hypothesis. Values such as -F(l, 4) and -F(2, 2) 
have no arithmetic interpretation, however they easily reduce to two-dimensional 
sums via classical theta series results. The usual method for detecting lacunarity, is 
to expand the associated cusp form in an infinite series. If one writes 



then the non- vanishing aj's should have zero-density. For the cases discussed herein 
it is usually necessary to compute thousands of coefficients to observe lacunarity. 
Additional non-obvious lacunary values include -F(2, 9) and -F(4, 7, 7, 28). It is nec- 
essary to employ eta function product identities to deal with these last two cases. 
By a result of Ramanujan [3, p. 210, Entry 56], we have 



Substituting classical theta expansions for ef, 62/61, and 6^/62 [10, pg. 114-117], 
leads to 



7l{q'')ri{q')ri{q')r]{q'') = q^-+^+^+^)l^^ (ao + a,q + a^q^ + . . . ) 



(14) 



3616269618 



00 



(2n+l)^ + (2fc + l)^ 
8 



Mq)r]{q^)r]{q')r}{q'') = - J](-l)"(2n + l)q 



(15) 




n=0 
k=0 




n=0 
k=—o 



00 



Because 616269618 is a finite linear combination of two-dimensional theta series, it 
must be a lacunary eta product. Formula (15) is the main ingredient needed to 
relate F[2, 9) to hypergeometric functions and Mahler measures. 
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Theorem 2. Let t = \fVl, then the following identity is true: 



144 



■F(2, 9) = - 3m (4z) + 2m (^-^ (4 - 2t - 2t^ + t^)^ 



(16) 257r2 

+ m(4i (7 + 4t + 2t2 + t^)) . 

Proof. The most difficuh portion of the calculation is to find a two-dimensional theta 
series for 616269618. This task has been accomplished via an eta function product 
identity. The remaining calculations parallel those carried out in [16]. Integrating 
(15) leads to 

3 ^. , ^, , (-l)"(2n + l) 



-F(2,9) + F(l,2) =4^ 

(17) 



25 ^ ' ' ■ ^ ' ' ^((2n + 1)2 + 9(2^ + 1) 

k=0 



_^ (-l)'^+'=(2n + l 



k=—oo 



There are two possible formulas for F(l,2) [15]: 

2 2 

(18) F(l, 2) =^m (2V2) = ^m (4^ . 

By the formula for -F(i,2)(3) in [16, Eq. 115], we also have 

(19) y _tl)!!!(^!l±l) =^^(4z(7 + 4t + 2t2 + t3)) 
^ ((2n + l)2 + 9(2A;)2)2 48 ^ ^ ' 

fc=— 00 

where t = \^12. Next we evaluate the remaining term in (17). Notice that for x > 

(-l)"(2n + l) 
^ ((2n + 1)2 -|-x(2A: + 1)2)2 

fc=0 

= '^j^u ||f^(-l)"(2n + l)6--("+^/2)'"j ||^g-^x(fc+i/2)2«j 
By the involution for the weight-3/2 theta function 

J](-l)"(2n + i)e-("+V2)^" = ^(-l)"(2n + l)e-^^^^^'^'^'"^ , 

n=0 ^ n=0 



this becomes 

oo 
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-l)"(2n+ 1) 



n=0 
k=0 



{{2n + ly + x{2k + lyy 



2 ^ roc 

n=0 
fc=0 

2 ^ -I \ 

^ V^z-^Nn lf^ljl_U g-^(2n+l){2fe+l) 



n=0 
fc=0 



* n=0 ^ 



X g-7rVx(n+l/2) 
g-7rv^(n+l/2) 



Applying formulas (1.6), (1.7), and (2.9) in [12], we have 



IT 

3V^ 



m 



m 4 



■i-/r 



a 



x/i 



1 + ,/l - a 



z/4 



where ax is the singular modulus (recall that Ux = 1 — V9^(— e ^^)/Lp'^[e '^^)). The 
second degree modular equation shows that 



and hence we obtain 



(20) 



1 - V 1 - Otx/A 

1 + a/1 - ax/4 
(-l)"(2ra + 1) 



TT 



^((2n + l)2 + x(2A; + l)2)^ 32Vi 

fc=0 



It is well known that a„ can be expressed in terms of class invariants if e Z: 

a„ = ^ (l - Vl - l/G^S^ 



and the values of G„ have been extensively tabulated, read: [4, p. 188]. Setting 
n = 9 yields 



V 



V3+ 1 



(1 - 4t + t^) 

(4 - 2t - 2t^ + t'f 
32 ' 
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where t = \fVl. It follows immediately that 

, , ^ (-l)"(2n + l) / 1 , 1 

21 V ^ '-^ ^ = — m — (4 - 2t - 2^2 + . 

^ ((2n + 1)2 + 9(2A;+ 1)2)2 gg J J 

k=0 

The proof of (16) can be completed by combining (17), (18), (19), and (21). □ 

In order to avoid tedious calculations we have chosen to omit the explicit formula 
for F(4, 7, 7, 28) from this paper^ It suffices to say that the sum reduces to an 
extremely complicated expression involving hypergeometric functions and Meijer 
G-functions. The key modular equation used to prove lacunarity is due to Somos 
[20, Entry 928,9,35]: 

(22) 28646^628 = -7eie? - + 8^614. 

By classical theta expansions [10, pg. 114-117], the identity can be rewritten as 

00 

28r7(g^)r/2(g7)^(g28) = - 7 ^ (-1)"+^X2A; + l)g" 

n=— 00 
k=0 



- J2 {Qn + l)q- 



n=— 00 
k=0 

00 



+ 8 ^ (-l)"+*^(3n + l)g- 

n,fc=— 00 

and as a result it is easy to see that 6467628 is lacunary. 

Apart from F(2,9), F(4, 7, 7, 28), and the examples discussed in [16], we are 
not aware of any additional lacunary values of F{a, b, c, d) (although they probably 
do exist). It is also possible to find two-dimensional reductions for certain linear 
combinations of lattice sums, however these formulas are generally less interesting 
than the previous examples. To give a single case let us briefiy consider the following 
modular equation [20, Entry X5o,6,8i]: 

e2 e2 

(23) 56162625650 + 26^62650 + 2616^625 = -6162 + el— + el — . 

625 650 

All three eta quotients on the right-hand side of (23) have two-dimensional theta 
series expansions. As a result it is possible to prove that 

^F(2, 25) + |f(1, 1, 2, 50) + (1, 2, 2, 25) 

(24) 



^ ' 5m {Ai) + 2m {4y/a^) + m ( a ' 



80 V V V V V "25 



""^Thc formula is available in previous versions of this paper on the Arxives 
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8 8 

where 0^25 = 2T3" {V5-1) (4^-1) . There are many additional results along the 
lines of (24) which we will not discuss here. 

3.2. Non-lacunary cases. In instances where F{a,b,c,d) does not reduce to a 
two-dimensional sum, the calculations become far more difficult. The recent proofs 
of formulas for -F(l, 5), F{2, 3) and -F(3, 5) are all based upon new types of g-integral 
transformations [17], [18]. The fundamental transformation used to prove a formula 
for F(2, 3) is 

Q 

where u = e^'^*/^. When x = 1 the left-hand side equals 4F(2, 3) (to see this use 
q^^^ijj{q) = yf{q^)/ri{q) and ^p{—q) = ff {q) / f]{(f)) ^ ^"^^ the right-hand side becomes 
an extremely complicated elementary integral. The most difficult portion of the 
calculation is to reduce the elementary integral to hypergeometric functions. It was 
proved with difficulty that 



/ill 1\ 

F(2,3) = -m(2) = -,f,(-yijj 



Boyd's numerical work was instrumental in the calculation, because it allowed for 
the final formula to be anticipated in advance. 

There are many cases where non-lacunary lattice sums reduce to elementary inte- 
grals, but where the integrals are extremely difficult to deal with. We recently used 
the method from [17] to find a formula for -F(l, 8): 
(25) 

zrn c^ 97r^/2 {I - kf + 2^2{k + k^) ^ /I + 2A; - P + 2v/F^\ , 
"^^'^'^ = ^1 (l + fc)(fc + fe3)3/4 '-^[ TTF 

We checked this monstrous identity to more than 100 decimal places by calculat- 
ing F(l,8) with (11). We can only speculate that the integral should reduce to 
something along the lines of (16). 

Occasionally eta function identities provide shortcuts for avoiding integrals like 

(25) . We have already demonstrated that linear dependencies exist between lattice 
sums (see (24)). In certain cases it is possible to relate new lattice sums to well- 
known examples. Consider a forty-fifth degree modular equation due to Somos [20, 
Entry 0:45,4,12]: 

(26) 6616569645 = -6165 - 263615 - 969645 + 63 + 5615. 

We were unable to prove (26) by elementary methods. Integrating (26) leads to a 
linear dependency between three lattice sums. We have 

(27) 9F(5,9) = 45F(1,1)-50F(1,5). 
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Both 1) and F(l,5) equal values of hypergeometric functions [15], [17]. Since 
the difficult task of evaluating F(l,5) is accomplished in [17], we easily obtain the 
following theorem. 

Theorem 3. Recall that n{a) is defined in (1). We have 

(28) = (^^) " 

There are various additional formulas which follow from Boyd's Mahler measure 
conjectures. A proof of Boyd's conductor 30 conjectures would lead to closed forms 
for both F{2, 15) and -F(2,5/3). To make this explicit we use two relations. First 
consider a four term modular equation which Somos highlighted in [19]: 

(29) 616365615 + 26266610630 = 6162615630 + 636566610. 

Integrating (29), and then using the evaluation F(3,5) = 47r^m(l)/15 from [18], 
leads to 

(30) F(2,15) + 4f(2,0 =^m(l). 

Next we require an unproven relation. Boyd conjectured^ that for a conductor 30 
elliptic curve 

L(i?30,2) = ^^7(3), 
15 

where g{a) is defined in (1). The modularity theorem guarantees that L(i?3o,2) = 
L^fsQ, 2), where /3o(6^''*'^) is a weight-two cusp form on ro(30). Somos has calculated 
a basis for the 1-dimensional space of cusp forms on ro(30), and consequently the 
cusp form associated with conductor 30 elliptic curves is given by 

/30(g) = v{q'Mq')v{qXq'')-v{qMq'Mq''Mq'')- 

Upon integrating this cusp form, Boyd's conjecture becomes 

(31) f(2,0-1f(2,15)^^^(3). 

Combining (30) and (31) leads to a pair of conjectural evaluations. 

Conjecture 4. Recall that m{a) andg{a) are defined in (1). The following formulas 
are numerically true: 

(32) i^F(2,i5) = 3m(l)-^?(3), 

(33) S^l'^'^^ =3m(l) + (7(3). 



7r2 



Tracking backwards shows that a solution of either (32) or (33) would settle Boyd's 
conductor 30 Mahler measure conjectures. Proofs remain out of reach, however we 
are optimistic that both identities may eventually be established using Eisenstein 
series identities contained in [1]. 



^See Tabic 2 in [6]. In our notation, Boyd's entries correspond to values of 17(2 — k) 
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4. Conclusion: Conductor 17 elliptic curves 

A strong connection exists between lattice sums and Mahler measures, however 
this relationship has limitations. While our ultimate goal is to "sum up" F{b, c) for 
arbitrary values of b and c, it is important to realize that this would only settle a 
small portion of the conjectures in Boyd's paper [6]. Conductor 17 curves are the 
first cases in Cremona's list [7], where L{E, 2) probably does not reduce to values of 
F{b, c). If we let Eu denote a conductor 17 curve (we used y'^ + xy + y = —x"^ — x), 
then 

17 

2^ 



(34) —L{E,r, 2) = m I ^ ) - m ( ^17 




We discovered (34) via numerical experiments involving elliptic dilogarithms. The 
cusp form associated with conductor 17 curves is stated in [9]. We have 

^ ^ ry(g2)^(gi7)^2(^68) r]{q)r]\q^)r]{q^^) ' 

Since L{En,2) = L{fn,2), formula (34) can be changed into a comphcated (!) 
elementary identity. There does not seem to be an easy way to relate L{En, 2) 
to Mahler measures of rational polynomials. We surmise that this is the reason 
conductor 17 curves never appear in Boyd's paper [6]. Given the complexity of 
fn{.(l)i we feel confident to conjecture that L{En,2) is linearly independent from 
values of F(6, c) over Q. 
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